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ON  PARAMETERS  OF  A  SUPERSONIC  COMBUSTION  CHAMBER  WITH  ORGANIZA¬ 
TION  OF  COMBUSTION  AT  THE  FLAME  FRONT 

E.  L.  Solokhin,  V.  A.  Mironenko,  V.  I.  Ivanov 

Notations  P,  T  are  static  pressure  and  temperature;  P» ,  T* 
are  stagnation  pressure  and  temperature;  M  is  the  velocity  coef¬ 
ficient;  J  is  the  total  flow  momentum  in  the  section  being 
examined;  F  is  the  area;  P.f  is  the  force  of  friction;  G  is  the 
consumption  by  weight;  k  is  the  adiabatic  index;  C„  is  the  ther¬ 
mal  capacity;  Q  is  the  heat  given  off  during  chemical  reaction; 

A  is  1/427-heat  equivalent  of  mechanical  work;  U  is  velocity; 
g  is  the  acceleration  of  gravity;  i«  is  the  total  enthalpy;  ot  is 
the  angle  of  inclination  of  the  flame  front  with  respect  to  the 
velocity  vector  of  the  fuel  mixture;  uT  is  the  turbulent  velocity 
of  flame  propagation;  p  is  density;  C,  is  the  drag  coefficient; 

O  is  the  coefficient  of  preservation  of  total  pressure;  £  Is  the 
friction  coefficient;  f7  is  the  perimeter;  9  is  the  width  of  the 
combustion  chamber;  Y*  is  the  halfheight  of  the  combustion 
chamber  at  intake;  is  the  ordinate  of  the  flame  front;  Y„  is 
the  current  halfheight  of  the  combustion  chamber;  Hu  is  the  lower 
calorific  value  of  the  fuel;  Ct  Is  the  fuel  weight  concentration; 
K*  is  a  coefficient;  a  is  the  velocity  of  sound;  0fM  is  the  excess 
air  ratio  at  the  Intake  in  the  combustion  chamber;  oc„  is  the 
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local  value  of  the  excess  air  ratio;  L0  is  the  stoichiometric 
mixture  ratio;  uH  is  the  normal  velocity  of  flame  propagation;  As 
(2.4)  is  an  experimental  coefficient;  U  is  the  pulsation  compo¬ 
nent  of  velocity;  P°  ,  T0  are  the  pressure  and  temperature  of  the 
mixture  for  the  prescribed  value  for  which  the  normal  velocity 
of  flame  propagation  is  equal  to  u0 H  (experimental  values);  £t  is 
the  degree  of  turbulence;  R  is  the  gas  constant;  Le  is  the  dis¬ 
tance  to  the  flame  line;  L0  is  the  length  of  the  flame;  x,  y  are 
coordinates;  A,  B,  B  are  typical  longitudinal  sections  of  the 
combustion  chamber;  "CM"  is  mixture;  "nM  denotes  reaction  pro¬ 
ducts;  **B"  intake  into  the  combustion  chamber. 

/ In  some  engineering  problems,  it  is  necessary  to  burn  fuel 
in  the  combustion  chamber  with  supersonic  flow.  As  a  rule,  the 
scheme  of  organization  of  the  process  in  such  a  chamber  presup¬ 
poses  a  separate  accompanying  feed  of  fuel  and  oxidant  in  which 
combustion  of  fuel  takes  place  in  a  diffusion  flame  front.  In 
this  article  we  give  theoretical  results  of  investigation  of  a 
supersonic  combustion  chamber  in  which  combustion  of  the  fuel 
mixture  takes  place  in  an  oblique  flame  front  stabilized  by  an 

external  source  (analogous  to  the  subsonic  cpmbustion  chambers  of 

Csvc.t* 

ramjets).  The  possibility  of  the  existence  of^an  oblique  flame 

"■  •'-'t  r  v  i  r.  ■  i 

front  in  a  supersonic  flow  of  fuel  mixture  was  proved  experi- 
menta Ilyin  C  5 ] . 

V 

Diagram  of  the  Combustion  Chamber  and  a  Model  of  the  Process 

A  diagram  of  the  combustion  chamber  is  shown  in  Fig.  1.  The 
following  organization  of  the  combustion  process  in  the  chamber 
is  proposed.  A  supersonic  flow  of  fuel-air  mixture  arrives  at  the 
combustion  chamber  intake  (section  A-A).  Fuel  gas  is  introduced 
into  the  same  section  A-A  which  has  a  static  temperature  higher 
than  the  static  temperature  of  the  main  flow.  The  fuel  gas,  on 


mixing  with  the  fuel-air  mixture,  ignites  the  latter  at  a  certain 
distance  from  section  A-A.  The  ignition  point  or  the  flame 
stabilization  line  (in  section  B-B)  is  adopted  as  a  constantly 
acting  ignition  source,  which  ignites  the  turbulent  oblique  flame 
front.  The  fuel  mixture,  on  passing  through  the  flame  front, 
undergoes  combustion.  The  method  of  obtaining  a  constantly  acting 
ignition  source  may  be  a  different  one,  for  instance,  a  pilot 
flame  in  the  bay  of  the  chamber  wall.  Conditions  for  obtaining  an 
ignition  source  and  conditions  for  it  to  be  stable  are  not 
examined  here. 


~ig.  i.  Combustion  chamber  diagrams  l  -  f  .s! -air 
ni  c.ure;  2  -  fuel  gas!  3  -  projection  of  the 
• ’ ~ ~e  stabilization  line;  4  -  reaction  products! 

5  -  f lame  front. 


A  System  of  Equations  for  the  Process  in  the  Combustion  Chamber 

In  putting  together  a  system  of  equations  describing  the 
process  in  a  supersonic  combustion  chamber,  the  following 
assumptions  are  made: 


1.  The  flame  front  is  infinitely  thin. 


2.  The  static  pressure  in  any  cross  section  of  the  chamber 
is  constant. 

3.  The  boundary  layer  of  the  chamber  wall  is  partially 
accounted  for  by  the  force  of  friction. 

4.  Dissociation  of  reaction  products  is  not  taken  into 
account. 

5.  The  field  of  parameters  of  the  gas  flow  in  a  cross 
section  at  the  chamber  Intake  are  uniform. 

The  adopted  assumptions  sharply  idealize  the  real  process  in 
the  combustion  chamber.  Nevertheless,  examination  of  the  problem 
in  this  formulation  can  give  aii  idea  as  to  the  order  of  the  major 
variables  characterizing  the  combustion  chamber  on  the  whole:  the 
necessary  length,  the  coefficient  of  preservation  of  full  pres¬ 
sure,  the  velocity  coefficient  and  other  parameters  on  the  outlet 
of  the  chamber. 

All  the  calculations  given  below  are  relative  to  a  combus¬ 
tion  chamber  of  rectangular  section  with  parallel  side  walls. 

The  problem  is  to  find  the  ordinate  of  the  flame  front  Y^s 
and  the  parameters  of  the  reaction  products  and  fuel-air  mixture 
along  the  length  of  the  combustion  chamber.  As  the  parameters  to 
be  determined  along  the  length  of  the  combustion  chamber  we  have 
chosen:  the  static  pressure  P,  which  for  the  reaction  products 
and  the  fuel  mixture  in  the  examined  section  is  the  same,  the 
velocity  coefficients  Mn  and  Me„,  the  velocities  of  sound  an  and 
ae n  .  The  quantities  are  sufficient  for  determining  any  other 
parameters  of  the  reaction  products  and  fuel  mixture. 


To  solve  the  problem  we  use  the  following  equations,  written 
out  for  an  arbitrary  section  of  the  chamber. 


1.  The  equation  for  preservation  of  total  momentum  taking 
into  account  friction  on  the  chamber  wall 


PdF-Pmp  . 

F, 


2.  The  consumption  equation 


G„  !  GlM  -  Gt 


3.  The  adiabatic  equation  for  the  fresh  fuel  mixture 


r.  \  r» ) 


4.  The  energy  equation 


Gti*t~G„l*u  |  GrMlmtM-~(}t 


5.  The  enthalpy  equation  for  the  fresh  fuel  mixture 


l\  -  1\m- 


6.  The  equation  determining  the  position  of  the  flame  front 


in  the  flow 


sin*  — 


Let  us  rewrite  equations  (l)-(6)  in  the  parameters  selected 
as  unknowns,  having  first  expressed  the  values  of  Uc „  and  uT  in 
terms  of  them. 


The  force  of  friction  on  the  segment  from  0  to  X  is  deter¬ 


mined  by  the  equality 


If  one  takes  into  account  that 


fl  =  «+2yfc  , 


(7) 


(8) 

(9) 


then  equation  (7)  can  be  put  in  the  form 


( 10) 


The  quantity  of  heat  Q  given  off  on  the  section  from  0  to  X 


»• 


Q  —  *//«j  pm  UlMC,  dy 

y 


or 


Q  =  *Huf<f»g{'^C1dy 
y*  a~ 


The  fuel  concentration  in  the  mixture 


C,  = 


•m  La  +1 


for 


and 


C,  = 


I 


Lo  +1 


for  1 


(11) 


( 12) 


( 13) 


( 14) 


Since  there  is  a  limited  amount  of  data  on  the  velocity  of 
flame  propagation  in  supersonic  flows,  we  use  the  equation  for  u 
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proposed  in  C 2 3 : 


.1. 


tr, 


«  + 


Uh 


k'K'+rf) 


( 15) 


where  the  normal  velocity  of  the  flame  propagation  uH  is  calcu¬ 
lated  by  the  formula  derived  according  to  the  data  of  [33, 


.  -U.J5 


“•=-  (tt)  (£r) 

In  equation  (15)  the  relative  Heating  up  &  is  equal  to 


( 16) 


(17) 


while  the  pulsation  component  of  the  mixture's  velocity  is 


.  (18) 

It  should  be  noted  that  the  choice  of  the  form  of  the 
equation  for  uT  in  this  statement  of  the  problem  is  not  really 
important.  The  important  thing  is  that  it  makes  it  possible  to 
take  into  account  the  Influence  on  uT  of  the  greatest  number  of 
physical  parameters  of  the  fuel  mixture  while  not  contradicting 
the  experimental  data  on  the  velocity  of  flame  propagation  ir 
supersonic  flow  available  in  the  most  up-to-date  literature.  For 
comparing  the  numerical  values  of  uT  obtained  from  (15)  with  the 
numerical  values  given  in  C43,  we  made  the  appropriate  calcula¬ 
tions.  The  results  are  shown  on  Fig.  2,  where  we  see  that  in  the 
region  of  lean  mixtures  for  turbulence  degree  IX  the  agreement  of 
the  calculated  and  experimental  values  of  uT  is  rather  good. 


Fig.  2.  Dependence  of  turbulent  veloc  ity  of 
flame  propagation  in  a  supersonic  flow. 

Key:  (a)  m/sec;  (b)  Calculated!  (c)  Experiment! 

id)  N/m. 

Other  equations  that  we  used  in  calculations  of  uT  (not  given 
here)  did  not  affect  the  final  results  qualitatively,  and  the 
quantitative  differences  were  minor*. 

So  far,  there  have  been  few  studies  made  to  clear  up  the 
level  of  intensity  of  turbulence  in  hi gh- temperature  supersonic 
flows.  In  the  calculations,  one  takes  £=0.01,  which  agrees  with 
the  data  derived  by  Baranova  in  experiments  with  cold  supersonic 
flows  with  M=2. 5  C5]. 

’Similar  results  are  obtained  if  one  takes  the  velocity  of  flame 
propagation  for  different  zones  of  the  flame  front  (the  forward, 
middle  and  trailing  zone  or  for  the  maximum-section  zone). 


One  has  gotten  roughly  the  same  level  of  intensity  of 
turbulence  in  experiments  carried  out  under  the  direction  of 
Tyul'panov.  Bearing  in  mind  that 

=  t^arcslna,  (19) 

dx 

in  final  form  equation  (6)  can  be  written  as  follows: 


dv,t 

,  dx 

—  i &  «r L^i|i  ■  ■■  ■ 

(ifM  Air. «f 

Vr.  A7CVx 

o  4  #r  &1'M  Mcm  (  U 

Y(poY0M  > 

'In  -  Kim  Ki  m 

(  ««u  v/v 

« A  #*  / 

Uim1 Ku  Kn 

]/  /„[ 

L(e^n 

(20) 


Accordingly,  equtions  (l)-<5)  as  functions  of  the  unknowns 
can  be  put  in  the  following  forms: 


«y« Pi (1+  Kcm  Af. 2 )  —  dP  Kf„  M\m) l-(y«  —  y»)X 

X  (H-/C,  M„  * )  ]  +  ~J  PMJ  (H-2 y*)dx-  jb  PdF 


>-B 


P  If  „  J.  D  f  l'<  »  Air  u  ,  .  Ki./VI,, 

p •  Km  a .  y-~  p  l *  1  (y *  i r  1- 


(21  ) 


(22) 


(23) 


'Pt  KcJ  - y*  Cp*  T.  # - \\cM'  gPMcjuUcM  A 

a§ 

J  )  y,P 

V.  Kr*  Rc m  2g  / 

CLC  “—  +  A-  M,,*  \y«  -->»■>)  - 

\K,  «n  / 


(24) 


(25) 


Cm  T« 


>,  &CMg  Rcm 


Thus,  to  determine  the  6  unknowns  (Y$,  P,  Mn,  Mc H  t  a„  ,  ac « ) 
along  the  length  of  the  combustion  chamber,  one  has  the  six  equa¬ 


tions  ( 20) - ( 25 ) . 


Results  of  Numerical  Solution  of  the  Derived  System  of  Equations 


The  system  of  equations  was  solved  on  a  digital  BESM-2M 
series  computer.  The  goal  of  the  solution  was  to  derive  the  dis¬ 
tribution  of  the  parameters  along  the  length  of  the  combustion 
chamber,  to  determine  the  length  and  position  of  the  flame  jet, 
to  clarify  the  special  features  of  the  flow  in  the  chamber 
depending  on  the  physical  flow  parameters  at  the  chamber  intake, 
the  chamber  geometry,  the  kind  of  fuel,  etc. 


The  chamber  chosen  for  study  is  symmetric  with  respect  to 
the  plane  passing  through  the  x-axis  perpendicular  to  the  y-axis 
(Fig.  1)  with  dimensions:  height  at  inlet  2y=0.05  m,  width  at 
inlet  ff=0.1  m,  the  law  of  variation  of  chamber  area  along  the 
length  is  given  by  the  equation 


F=2*(y.  +KPX) 


(26) 


Some  obtained  results  are  shown  on  Figs.  3-8,  from  which  the 
following  inferences  may  be  drawn: 


1.  To  obtain  short  flame  Jets,  it  is  necessary  to  use  fuels 
having  high  velocities  of  flame  propagation  (for  example, 
hydrogen  +  air).  Ue  see  from  Fig.  7  that  the  flame  jet  length  for 
a  kerosene-air  mixture  is  4-5  times  longer  than  for  a  hydrogen- 


air  mixture. 


2.  The  flame  jet  length,  and  consequently  also  the  length  of 
the  combustion  chamber  for  broad  ranges  of  variation  of  the 
physical  parameters  at  the  chamber  inlet,  are  not  long  and  of  the 
same  order  as  for  subsonic  combustion  chambers  (see  Figs.  3,  4, 

5,  6,  on  which  the  dependence  are  given  of  the  flame  jet  length 
on  the  number  M,  the  stagnation  temperature,  the  excess  air 
ratio,  the  degree  of  expansion  of  the  chamber  and  the  inlet 
pressure  ino  the  combustion  chamber). 

3.  To  burn  the  greatest  amount  of  high-calory  fuel  (for 
without  choking  the  chamber,  it  is  necessary  to  make  it 

expansible  (see  Fig.  6).  And  the  degree  of  expansibility  (which 
corresponds  to  the  coefficient  Kr )  must  be  the  greater  the 
smaller  the  number  M  at  the  chamber  intake.  From  Fig.  6  we  also 
see  that  the  degree  of  combustion  chamber  expansion  barely 
influences  the  flame  jet  length. 

4.  The  greatest  Influence  on  the  flame  jet  length  turns  out 
to  be  the  excess  air  ratio  a  (see  Fig.  5).  The  stagnation  tem¬ 
perature  T,  antj  the  static  pressure  p  at  the  chamber  intake 
are  of  less  influence  than  <x  (see  Figs.  3,  4). 

5.  Calculation^  have  shown  that  for  appropriate  combinations 
of  the  flow  parameters  at  the  chamber  intake,  the  chamber 
geometry  and  the  kind  of  fuel,  one  of  the  following  two  typical 
flow  regimes  results: 

a)  AI„  >  1,  Mcm  >  1. 

b)  Af„<l, 

l.e.,  in  a  section  of  the  combustion  chamber  there  can  simulta¬ 
neously  exist  both  supersonic  and  subsonic  regions  of  flow 


(Fig.  8).  It  is  characteristic  that  the  switchover  from  one 
regime  of  flow  to  another  under  variation  of  any  parameter  is 
sa 1 t i f orm. 


In  conclusion,  note  that  the  presented  theory  of  the  process 
in  a  supersonic  combustion  chamber  and  the  numerical  results  need 
to  be  checked  experimentally  and  sharpened,  since  not  very  much 
is  now  known  about  the  values  of  some  of  the  basic  parameters 
<£t,  ur )  for  supersonic  flows. 
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Fig.  3.  Dependence  of  flame  jet  length  on  :ne 
number  M  and  pressure  P  .  Mixture  Ha  +  air, 
TV=2000«K,  6T=0.1,  K*  =0.3. 

Key:  (a)  N/ma . 
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Fig.  4.  Dependence  of  flame  Jet  length  on 
temperature  T,  •  and  number  Me  .  Mixture  H2  +  air, 
*«  =  1.  P«,  =2943  N/ma  ,  ST=0.01,  K,=0.3 {  1  -  region 
of  stall-free  flow;  2  -  region  of  choking. 
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Fig.  5.  Dependence  of  flame  jet  length  on 
0Lff  and  M*.  Mixture  H,  +  air,  TV=2000*>K. 

P$  =2943  N/m*  ,  a.«-0.01,  Kf  =0.2,  1  -  region  of 
stall-free  flow;  2  -  region  of  choking. 


•i  < 


yf‘ 


'•■KV; 


t’*  r- 


Fig.  6.  Dependence  of  flame  jet  length  on 
K*  --nd  M  .  The  mixture  H*  +  air,  Ts*=2500“K 
P*=4905  N/ma ,  a.*=l,  £T=0.0i,  i  -  region  of 
stall-free  flow;  2  -  region  of  choking. 
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Fig.  7.  Dependence  of  flame  jet  length  cn  M 
for  different  fuel  mixtures,  i  -  kerosene  +  air 
2  -  Ha  +  air,  P  =4905  N/m*  ,  T  =  2500«K,  Ql*  =  i, 

K*  =0.2,  £T  =0.  01 . 
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Fig.  8.  Dependence  of  Mrt,  Mc „  and  P  on  K„ 
at  0.3  m  from  the  intake.  Kerosene  +  air,  M* 
P*  =4905  N/m*  ,  Tj  *  *2000*  K ,  atf*i.4,  £■,  -  9i .  e  <B1 . 
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